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|V] ; ABSTRACT 

. The time-scale over which and the modahty by which young stehar objects (YSOs) 

disperse their circumstellar discs dramatically influences the eventual formation and 
evolution of planetary systems. By means of extensive radiative transfer (RT) mod- 
elling, we have developed a new set of diagnostic diagrams in the infrared colour-colour 
O ■ plane (K- [24] vs. K- [8]), to aid with the classification of the evolutionary stage of 

YSOs from photometric observations. Our diagrams allow the differentiation of sources 
c/3 ■ with unevolved (primordial) discs from those evolving according to different clearing 

^ ■ scenarios (e.g. homologous depletion vs. inside-out dispersal), as well as from sources 

that have already lost their disc. Classification of over 1500 sources in 15 nearby star- 
forming regions reveals that approximately 39 % of the sources lie in the primordial 
J> \ disc region, whereas between 31 % and 32 % disperse from the inside-out and up to 

OO . 22 % of the sources have already lost their disc. Less than 2 % of the objects in our 

^>0 ■ sample lie in the homogeneous draining regime. Time-scales for the transition phase 

I are estimated to be typically a few 10^ years independent of stellar mass. Therefore, 

regardless of spectral type, we conclude that currently available infrared photometric 

■ surveys point to fast (of order 10 % of the global disc lifetime) inside-out clearing as 

the preferred mode of disc dispersal. 
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1 INTRODUCTION iMuzerolle et all (|2010l ) by u sing IRAC data a n d the 24 ^m- 

„ , , , ^ , fluxes obtained with MIPS (iRieke et al.ll20o3 '). lFurlan et all 

Protoplanetary discs are an unavoidable consequence or the innr,A\- .-, ^ .^ ^■ i -,1,1 

(I2009l )investigated the disc evolution m the three promi- 

star formation process and harbour the material from which 



planets may form. For this reason, it is essential to under- 
stand their properties, evolution and final dispersal. For a 



nent star-forming regions Taurus, Ophiuchus, and Cha I 
and found similar transitional disc fractions in all three re- 



gions. This study is based on spectro scopic data taken w ith 
^eceM^^glgtggetary discs and their evolution ^^^^^^^^ Spectrograph' jHouck et al.ll2004h . A 

— ; ; ~ — spectroscopic study of the young s tellar objects in Taurus 

The detection and characterisation or circumstellar , , 1 ,• , ,1 [7; — , — ; — Tl JT^TTTHn ■^• r ^ 

has been published by llurian et al.l (120111 ). Ine mid-mrrared 
discs was a primary goal or the Spitzer Space Telecope . , „ / -^^ 1 i- 1 ■ j_- ^ 1 -^1 

I 'ni 1'^ ° f f f variability o r (pre-)transitionai di sks was investigated with 

^Werner et alj | 2004^ ) and successfully featured in the Legacy ^^.^^^^ re.g. lEspaiUat et all j20Tlh ). 
Prog ramme From Molec ul ar Lores to Planet -Formmg Dust 

fc2d: lEvans et al.l (|2009h '). lAUen et all (|2004l ') compared the The interpretation of infrared colours in relation to 

predicted colours of young stellar objects with the photomet- the evolutionary state of a disc is, however, far from being 
ric data obtained w ith the 'Infrared Array Camera' (IRAC; trivial. This is particularly true with regards to the clas- 
iFazio et all (|2004| ')V Transitional discs in different young sification of transition discs, the latter being intended as 
clusters and associations were identified and classified by objects caught in the act of dispersal, i. e. going from disc- 
bearing to disc-less stars. The modality and time-scales of 
disc dispersal are especially relevant to the planet forma- 
* koepferl@usm.lmu.de tion process. In particular, inside-out dispersal from internal 
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photoe vap oration (e.g. Clarke et al. | 20pi| ): Alexander et al 
2006al lbl): lOwen et all (|2010l . l2011al lbl. |2012| ): ICorti et al 



20091 )^ could also provide a mechanism to stop migration 



and effe ctively 'park' planets at a par ticular location in a 
disc fe.g. lAlexander fc Armitagd (|2009h '). 

The evolution of the dust component in a disc 
is mirrored by the evolution of colours in the in- 
frared plane. By means of radiative tra nsfer mo delling, 
lErcolano. Clarke, fc Hall (2011, henceforth TeCH 111 ) identi- 
fied the regions in the K - [8] vs. K - [24] plane where pri- 
mordial discs, discs with inner-holes (i. e. presumably being 
dispersed from the inside-out) and discs which lose mass 
homogeneously at all radii, are expected to be found. Their 
study, which was limited to M-stars, showed that in the case 
of the cluster IC348, most discs disperse from the inside-out 
and undergo the transition on a short time-scale, as pre- 
dicted by standard photoionisation models. 

These conclusions a re in contrast with the conclusions of 
ICurrie fc Ken"^ (I2OO9I ). who claimed instead a large num- 
ber of 'homogeneously depleting' discs, for the same cluster. 
Such discrepancies highlight the need for detailed modelling 
in the interpretati on of IR colours of discs. 

The study of lECHllI was restricted to M-stars in only 
one cluster, which prevented the authors from being able to 
make a more general statement with regards to disc disper- 
sal. 

In this study we significantly improve on the work of 
IECHIII by performing further RT calculation to evaluate 
evolutionary tracks in the IR colour plane for stars of differ- 
ent spectral types. We then apply our results to the photo- 
metric data of 15 nearby star-forming regions, that we col- 
lected from the literature, in order to address the question 
of what is the preferred mode of disc dispersal. 

In Section [2] we describe our methods and the results 
from our radiative transfer calculations and present new evo- 
lutionary diagnostic diagrams. In Section [3] the collected ob- 
servational data are briefly described and analysed using the 
diagrams presented in the previous section. A short sum- 
mary and conclusions are given in Section 2) 



2 THEORETICAL INVESTIGATION 

Emission in a given wavelength band is produced by dust 
grains of given temperatures and sizes. The temperature of 
grains in a disc for a given central star temperature is de- 
termined roughly by the grain location within the disc and 
by the intervening material. Changes in the distribution of 
material within a disc are therefore reflected in its emer- 
gent continuous spectral energy distribution (SED). On this 
basis, the evolution of infrared colours of an object can be 
used to trace the evolutio n of its dust density distribution. 
We extended the work of lECHl l| and calculate theoretical 
evolutionary tracks in the K - [8] vs. K - [24] plane by means 
of radiative transfer calculations for discs of different geom- 
etry and different evolutionary stage surrounding stars of 
different spectral types. 



2.1 Modes of Disc Dispersal 

We consider three modes of disc dispersal: (i) homoge- 
neous draining (HD), (ii) inside-out clearing (IOC) and (iii) 



initial disc 



disc-less 




IOC 



K-[24] 

Figure 1. Expected tracks of discs depleting via the three dif- 
ferent modes of disc dispersal explained in Section l2.1l The OIC 
track is marked with a question-mark, since it has been drawn ac- 
cording to the expectations from the simple thought experiment 
described in Section [2. II But as will be discussed in Section 12.3.11 
the real tracks behave somewhat differently. 



outside- in clearing (OIC). The expected tracks in the in- 
frared two-colour plane are sketched in Figure [T] 



2.1.1 Homogeneous draining 

HD refers to a disc whose surface density decreases homo- 
geneously at all radii as a function of time, such that the 
surface density distribution E(t) is self-similar at all times. 
Such a mode of dispersal would be expected if viscous drain- 
ing were to act alone, undisturbed until all of the disc mate- 
rial is accreted onto the central object. In this case the total 
mass of the disc evolv es roughly as Mdisc oc E(t) oc t~^'^ 
iHartmann et allll998l ). In this scenario one would expect 
discs to evolve diagonally on the K - [8] vs. K - [24] plane 
(see Figure [ij as both the 8 /im and the 24 ^m band fluxes 
decrease roughly at the same rate as material is removed 
homogeneously at all radii. 



2.1.2 Inside-out clearing 

IOC refers to discs that disperse via the formation of an 
inner-hole (classical transition discs), which becomes larger 
and larger until all material is removed. There is currently 
a healthy debate in the literature with regards to possi- 
ble mechanisms that may be responsible for the formation 
of discs with inner-holes, of which a large variety is ob- 
served. Some of these mechani sms include pho toevapora- 
tion by the central obje ct fe.g. lAlexander et al.l (2006 a g); 
I O well et al] (|2010l . l2011al )), planet formation Ce.g.rZhu et"al] 
(|201ll ')'). grain growth Ce.g.lDuUemond fc Dominikl (|2005l )T 
close-binary effects fe.g. llreland fc KrausI ( 20081 )') and exter- 
nal photoevaporation by nearby high-mass stars (e.g. Scally 
fc Clarke 2001). The picture that seems to be currently 
emerging is that of multiple populations of inner-hole ob- 
jects, some of which may o r may not represent objects in 
transition (e.g. lOwen et al.l (|2012l i: lOwen fc Clarkel I2OI2I ')). 
It is beyond the scope of this work to explore these differ- 
ent possibilities; here we will limit ourselves to calculating 
SEDs of discs, which, by whatever mechanism, develop in- 
creasingly larger inner-holes. The IOC tracks are expected 
to proceed in a left-facing 'L'-shape (see Figure [1]) as the 
warmer material in the inner disc is removed first, causing a 
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Table 1. Spectral type intervals and stellar properties. 



index 


SpT 


T4K] 


R4Rq] 




3013427 


K4 


4590 


1.35 


1.21 


3014686 


K8.5 


4060 


3.73 


0.75 


3002746 


M4.5 


3305 


0.844 


0.228 


3014449 


MS 


2805 


1.8 


0.1 



Table 2 . Input parameter ranges use d in our simulations 
(see also lWhitnev et al.l l l2003al lbh and lRobitaille et al.l i200dt) 
for a detailed description). 



drop in the 8 fim flux, and the cooler material, responsible 
for most of the 24 fj,m emission, is removed afterwards. 



2.1.3 Outside-in clearing 

OIC refers to discs th at are eroded from the outside-in by 
a com p anion star (e . g.lBouwman et al.l (|2006l ): iMonin et al.l 
|2003); ICieza et all (|2009l ')') or by radiation from exter- 
nal sources. External photo-evaporation is for example ob- 
served i n the case of prop l yds in the Orion Nebular Clus - 
ter (e.g.lRicci et al.1 (l2008f ): iMcCaughrean fc O'Delil (|l99d ): 
IScallv fc Clarke! (|200ll ))." Whether the entire mass of a disc 
can be dispersed via this mechanism remains, however, very 
uncertain. It is also uncertain what percentage of discs in 
OB associations and clusters are exposed to strong enough 
radiation from nearby massive stars for this process to be 
significant. In any case, it is an interesting exercise here to 
explore the evolutionary tracks for such a population on the 
infrared colour-colour plane. One would naively expect an 
upside-down 'L'-shape this time, almost a mirror-image of 
the IOC tracks (see Figure[T|, as the 24 /^m flux should drop 
before the 8 yLxn flux. This is, however, not necessarily the 
case, as will be discussed in Section [2.3. II 



2.2 Numerical Set-up 

We use the 2 -dimensional M o nte Ca rlo dust radiative trans- 
fer code of IWhitnev et all (|2003al lbl). which is routinely 
applied to the model ling o f the SEDs of proto planetary 
discs (e. g. [RobitaiUc et al.l (|2006l , |2007| ), lECHllI ). We de- 
rive evolutionary tracks for stars of spectral types K4, 
K8.5, M4.5 and M8, using the stellar atmosphere m odels 
l|Castelli fc Kuruczl ll200i):lBrott fc Hauschildtl (|2005M from 
the librarv used bv lRobitaille et aL ( 20061 ) and available on- 
line at http://caravaii.astro.wisc.edu/protostars/. Ta- 
ble[T]summarises the stellar parameters and provides indices 
for the atmosphere files used. Table [2] lists the other input 
parameters used in the 1287 radiative transfer models. A 
detailed description of the geometrical set-up and the input 
parameters is given in iRobitaille et al.l (|2006l l. 

We follow IECHIII and consider several disc geometries 
exploring flrst the evolution of razor-thin discs for difi'erent 
dispersal modes and then extend our study to discs of finite 
thickness (Tabled. 



2.3 Results 

2.3.1 Razor-thin discs 

Figure [5] summarises the tracks for razor-thin discs in the 
four spectral type intervals considered here. Point F are face- 
on (cosi = 0.95) and point E are edge-on (cosi = 0.05) disc 



parameter 


range 


sequence 




10"^ M^, 


initial, IOC, OIC 




(10-2 - 10-«)M* 


HD 


rdiaci'max) 


100 AU 


initial, IOC, HD 


Tciiscimax) 


(0.02 - 100) AU 


OIC 






initial, OIC, HD 




1—30 Rsub 


lOC^ 






IOC 3 




0.0001 - 1.0 


initial, IOC, OIC, HD 




2.0 - 2.25 


initial, IOC, OIC, HD 


Bl 


1.0 - 1.25 


initial, IOC, OIC, HD 




10-10 


initial, IOC, OIC, HD 


Mass infall rate_^ 


lO-l^Mo/yr 


initial, IOC, OIC, HD 


Aperture 


100000 AU 


initial, IOC, OIC, HD 



" sublimation radius 
^ razor-thin 
finite thickness 

disc scale-height, see also Table |3] 
^ disc density exponent (~ r~°') 

f disc scale height exponent (~ r''); see also Table [3] 
9 cfli'cctivcly no viscous heating 
^ effectively no envelope 



inclinations. The points between the F and E represent in- 
termediate disc inclinations in steps of A cosi = 0.10. Point 
P shows the approximate location of disc-less objects. The 
stars mark the HD tracks for face-on and edge-on inclina- 
tions (green and blue lines, respectively). The triangles mark 
the IOC tracks for face-on and edge-on inclinations (yellow 
and cyan lines, respectively). The diamonds mark the OIC 
tracks for face-on and edge-on inclinations (black and purple 
lines, respectively). 

As expected (discussed in Section [2.1.11 and sketched in 
Figure [T|), the HD tracks follow a diagonal trajectory from 
points F and E to P, while the IOC tra cks show the typi- 
cal left-facing 'L'-shape also reported bv lECHlH for M-star 
models. The OIC tracks in Figure (2] however, do not show 
the upside-down 'L'-shape suggested by our simple-minded 
thought experiment (see Section r2.1.3|) . This is because the 
24 fim emission is not dominated by material at large radii 
in the disc, but contains a significant contribution from ma- 
terial in the inner disc that is located deeper inside and is 
hence cooler. 



2.3.2 Discs with finite-thickness 

We follow lECHllI and set-up discs of finite-thickness and 
nine dif i'erent geometries; we extend, however, the work of 
lECHllI to all the spectral types considered here. The re- 
sulting infrared colours are plotted in Figure |3] for all cases 
where disc extinction Av < 10 mag. The parameter set-up 
is summarised in Table |3] 

We explore three different cases of dust settling (by 
varying the parameter Hd) and three different levels of fiar- 
ing, in order to cover the parameter space occupied by pri- 
mordial discs of different geometries. Discs with Hd ~ 0.1 
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Table 3. Modelling parameters of primordial discs. 

flaring mixing 



disc kind 


B 


a 


b 
-— s> 


c 


Hd 


d 


e 


f_ 


razor-thin 


1.0 


X 






0.0001 


X 






ultra-settled 


1.0 


X 






0.1 


X 






optically thick 


1.0 
1.0 


X 
X 






0.5 
1.0 




X 


X 


ultra-settled 


1.13 




X 




0.1 


X 






optically thick 


1.13 
1.13 




X 
X 




0.5 
1.0 




X 


X 


ultra-settled 


1.25 






X 


0.1 


X 






optically thick 


1.25 
1.25 






X 
X 


0.5 
1.0 




X 


X 



constant opening i 

^ flared 

strongly flared 
strongly settled 

^ settled 

^ mixed 



are named 'primordial ultra-settled', while those with a 
larger Hd ^ 0.5 are named 'primordial optically-thick' discs. 
We note that Hd also influences the inner wall height, which 
may be the main source of variation of the NIR/ MIR 
colours. 

By inspection of the locus occupied by the evolutionary 
tracks of primordial (ultra-settled and optically-thick) discs, 
we deflne boundary boxes in the two-colour plane. Our re- 
sults are summarised in Figure|31 As for the razor-thin case, 
we carried out radiative transfer calculations to track the 
colour evolution of each of our primordial (ultra-settled or 
optically-thick) discs as they disperse via HD and IOC, al- 
lowing us to further identify the expected locus of objects in 
such evolutionary stages. For the sake of clarity, the tracks 
are not plotted in Figure [31 where only the derived bound- 
ary boxes are shown. In particular the boxes are labelled as 
A, B, C, D, E for 'primordial optically-thick' discs, disc-less 
sources, 'primordial ultra-settled' discs, IOC and HD respec- 
tively. Due to the overlap in the evolving discs regimes D and 
E, we cannot discriminate between HD and IOC for those 
objects observed in DE. As will be shown later, however, 
only a very small minority of the observational sample lies 
in this regime. 

Additionally, we run a set of models to allow us to con- 
strain the locus of envelope objects and estimate the con- 
tamination of the IOC objects by this sample . We se t up 
576 objects with envelopes (see IWhitnev et al] (|2003al )) for 
a selection of cavity angles (10° to 90°) and envelope infall 
rates Mq / yr to Mq /yr). Objects with envelopes 

and high infall rates show the greatest overlap with the IOC 
region, whereas objects with less infall are bluer. We have 
indicated the locus of these models with a dashed line in Fig- 
ure [3] above which envelope objects lie. By comparing with 
observational data (see Section [3} , we found, that primor- 
dial discs with envelopes might introduce a contamination 
in the IOC statistics of at most 10%. This bias has to be 
kept in mind when considering the raw statistics of the up- 
coming section. The boundary points of each box are listed 
in Appendix 1X1 



3 ANALYSIS AND DISCUSSION 

Our classification scheme, based on the radiative transfer 
modelling described in Section 12.21 allows us to use cur- 
rently available infrared photometric surveys from nearby 
star-forming regions in order to address the question of what 
is the preferential mode of disc dispersal. We have collected 
K-band, 8 fira and 24 /xm band data, as well as spectral types 
and extinction measurements of YSOs in 15 nearby star- 
forming regions. Table 31 lists the references for the data and 
summarises some of the physical properties of the regions. 

We applied our colour-colour diagnostic diagrams to 
classify 1529 objects in 15 nearby star-forming regions. The 
number statistics in each evolutionary stage and for each 
spectral type regime are presented in Table [5] It is worth 
noticing here, that we deal in most cases with small num- 
ber statistics; for this reason we provide actual numbers in 
Table [5] rather than percentages. In what follows, however, 
we prefer to use percentages in order to simplify the dis- 
cussion. The standard statistical errors (e.g. \/N), however, 
should be kept in mind when considering the raw percent- 
ages given. Given the overlap between regions D (IOC) and 
E (HD), discussed in the previous section, we can only pro- 
vide upper and lower limits to the count sources in these 
regimes. 

In summary, 38.9 % of the objects out of the entire sam- 
ple lie in the primordial disc region whereas between 31.3 % 
and 31.9 % disperse their discs from the inside-out and up to 
21.7% of the sources have already lost their disc. So, almost 
a third of the available sources are currently clearing their 
discs from the inside-out. Only between 1.5% and 2.0% of 
the objects lie in the homogeneous draining region E. We in- 
terpret this result as strong evidence against homogeneous 
disc depletion as the main disc dispersal mode. 

Figures |4l to [I9l show the results for individual star- 
forming regions, which are also briefly summarised in the 
next section. The regions are sorted by age. 



3.1 Analysis 



Taurus (Figure Bl [5l) 
Myr (ILr 



For this 1 Myr (jLuhmanl l2004l : iLuhman et al] l2010l : 
iReipurth. B]|2008al) old reg ion w e have collect e d dat a both 
from iLuhman et al I (|201oi ) and iRebuU et al] (|2010l '). The 
statistics differ when using the two compilations as detailed 
in what follows. 51.2% and 41.6% of the Luhman and the 
RebuU discs, respectively, are still primordial optically- 
thick. None of the Luhman objects are consistent with being 
primordial ultra-settled discs, while 0.8 % of the RebuU 
sources fall in this category. 23.2% and 31.5% of objects 
in the Luhman and RebuU compilations, respectively, have 
already lost their discs. Currently, 13.7 % of the Luhman 
sources are depleting their discs via the IOC, while between 
17.9% and 19.1% of the RebuU objects are undergoing 
that process. Only one disc (0.5 %) of the Luhman objects 
and between 2.3% and 3.5% of the RebuU sources lie in 
HD region. In the RebuU case we also considered upper 
and lower limits in the 24 jim band as well as candidates by 
iRebuU et all (|2010l ). 

• NGC206 8/71 (Figure H 

In the 2 Myr (jFlahertv fc MuzeroUell2008l ') old star-forming 
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Table 4. Overview of the used observational data. 



region 


Ntot 


age [Myr] 


distance [pc] 


JHK 


[8] 


[24] 


Ax 


SpT 


Taurus^ 


211 


ll,2,3 


140* 


5 


S 




6 


2,6 


Taurus^ 


257 




137' 


5 


8 


8 


8 


8 


NGC2068/71 


27 


2" 


400'-' 


a 


!) 


a 


a 


a 


NGC2264 


313 


2l0 


760-913^^'^^ 


13 


13 


13 


13 


13 


Chal 


29 




162^'' 


1« 


1« 


IB 


17 


17 


IC348 


240 






•21 


21 


21 


21 


21 


NGC1333 


39 


<3^^ 


260^*^ 


m,25 


ili,25 


i«,25 




25 


Tr37 


57 


427,28,29 


900^" 


is 


2d 


iO 


28 


28 


Serpens 


69 


2-6^^ 


260-415^^'^^ 


34 


34 


34 


31 


31 


LupIII 


82 


2_g3b,3b 


200"'^ 


34 


34 


34 


36 


36 


Chall 


27 


4_53V 


178^'' 


34 


34 


34 


37 


37 


NGC2362 


10 


538 


1480" 


14 


14 


14 


39 


3a 


OBlbf 


13 


540 


440-41 


5 


42 


42 




43 


UppcrSco 


133 


544 


145**^ 


5 


4a 


47 


47 


47 


T] Cha 


14 




1004« 


5 


2 


2 


1 


4a 


250ri 


8 


7-10''" 


330='^ 


5 


i2 


42 




43 



L: Data from[Luhman ct al. (2010) 



R: Data from iRebull et al., (i201Q ) 



LuhmanI ||2004|) 

Luhman et al. I J 2OIOI) 



Reipurth . B. (2008a 
Kenvon ct al. ( 1994 
Cutri et al. (2003 



Luhman et alj (|2009|) 



Torres et al. I 2007 ) 
RebuII et atl 1I2OI0I) 



Flaherty fc MuzeroUj 1 I2OO8 
Teixeiral ||200S|'| 



Sung et al.l lll997l'l 
Baxter et al.l j2009l\ 



Teixeira et alj ll2012h 
Currie et al.l ^2009^ 



Whittet et al I lll997h 



Gutermuth et al .1 l|2009l ') 
LuhmanI (120071 ') 



Herbid il99^) 

Muench et'aU ( |2003|) 
Luhman et al.l (|20(33|'| 
Currie fc Kenvonl ( l200g|) 
HerbstI (I2OO8II 
Scholz ct alT (|l999l 'l 

Cernis ( 1993) 

Winston et al] (|2009|) 
Hirota et alj (I2OO8I) 



Sicilia-Aguilar et al 



Sicilia-Aguilar et al 



Sicilia-Aguilar et al 



2004 ) 



Contreras et alj (l2002l i 



Oliveira et al 



Straizvs et al 



200E) 



^ 20061) 



Dzib et alJ J2OI0I) 
Evans et al.l ( |2009^ 



Reipurth. b] (l2008bl) 
Men'n et air ( l2008l) 
Spezzi et all (I2OO8I) 



Maync & Naylor' ('2008') 
Dahm & Hillenbrand (20^) 
Hernandez et al. ( 20061) 
Hernandez ct al., (^20051) 
Hernand ez et aO (l2007bl ) 
Hernandez! (I2OIII) 



Hernandez et al. 



de Zeeuw et al.l jl999l)" 



(l2007al) 



Carpenter et al 



Carpenter et al 



( 2006 ) 
J2OO9I ) 



Mamaiek et 81.1(119991) 



Luhman fc Steeghj (|2003 ) 



Briceho et al. 
Briceho et al.l 



2007 ) 



20051) 



region NGC2068/71 81.5 % of the discs are still classified as 
primordial optically-thick. 14.8 % of the sources are clearing 
their discs from the inside-out. 

• NGC2264 (Figure [Tj 

25.2% of the discs of the 2 Myr (|Teixeira et al.ll2012l ) old 
cluster NGC2264 are still primordial optically-thick. 0.6 % 
are consistent with them being primordial ultra-settled 
discs. 9.9 % have dispersed already and between 56.2 % 
and 57.2 % are clearing via the IOC process. Only 0.3 % to 
1.3% of the discs are found at the edge of the HD regime. 
We include in this statistics also objects with only upper 
limit measurements in the 24 fim band. All the objects 
are classified as ca ndidates by radial velocity studies by 
iFtiresz et all (|2006l ). 

• Chal (Figure© 

In Cha I, a 2-3 My r (|Luhmanl l2007l : iReipurth. b1 l2008bl : 
iLuhman et aI.ll2O10l ) old region, 75.9 % of the discs are still 
primordial optically-thick. One source (3.4%) is clearing its 
disc from the inside-out, while another disc (3.4%) can be 
found at the edge of the HD regime. 



• IC348 (Figure [9| 

23.8% of the objects of the 2-3 M yr (|Herbig| Il998l : 
I Muench" et al.ll2003l : lLuhman et~aLll2003l ) old cluster IC348 
are still primordial optically-thick. 0.8 % are consistent with 
being primordial ultra-settled discs. 3.8 % have dispersed 
already and between 69.6 % to 70.0 % are clearing via 
the IOC process. We include in these statistics also ob- 
jects with only upper limit measurements in the 24 fim band. 

• NGC1333 (Figure [ini) 

66.7 % of the discs in the less than 3 Myr (|Winston et al.l 
l2009h old cluster NGC1333 are still primordial optically- 
thick. 28.2 % of sources are depleting their discs via an IOC 
process. 

• Tr 31 (Figure [TT1) 

In the 4 Myr (| Sicilia-Aguilar et al.l |2004 |2005| . |2006| ) 

old star-forming region Tr 37 89.5 % of the discs are 
still primordial optically-thick. None of the objects have 
primordial ultra-settled discs or have depleted their discs 
already. 7.0 % are currently clearing via the IOC pro- 
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cess. For Tr 37 no sources are found in the HD region. 
Note, that we also included YSO candidates in our statistics. 

• Serpens (Fi gure | 12ll 

For the 2-6 Myr (|01iveira et al.|[2009l ) old cluster Serpens we 
have classified 40.6 % of the discs as primordial optically- 
thick and one disc (f.4%) as primordial ultra-settled. 
18.8 % of the objects have already lost their discs. Between 
seven (10.1 %) and eight (11.6 %) discs are undergoing IOC. 
Between seven (10.1%) and eight (11.6%) discs lie in the 
HD regime. Candidates are included in the statistics. 37.5 % 
of the objects i n the HD box are classified as candidates by 
lOliveira et all (|2009l ). 

• Lup III (Figure [T3)l 

39.0 % of the disc s in the 2-6 Myr l|Reipurth. B.I l2008bl : 
iMerfn et all 120081 ) old star-forming cluster Lup III are 
still primordial optically-thick. One disc (1.2%) has been 
classified as primordial ultra-settled. 14.6 % of discs are 
undergoing IOC, while 31.7% have already dispersed their 
discs. We also included YSO candidates and sources for 
which only upper limits in the 24 fim flux are available in 
our statistics. In Lup III 7.3 % of the discs are found at the 
photosphere edge of the H D regime. These so urces are still 
classified as candidates bv lMeriri et al.l (|2008D . 

• Cha II (Fig ure [Till 

In the 4 -5 Myr (jSpezzi et al.ll200^ old star-forming region 
Cha II 63 % of the discs have been classified as primordial 
optically-thick, while 25.9 % have already depleted their 
discs. 7.4 % of the objects are clearing their disc from the 
inside-out. For Chall no discs are found in the HD region. 
We included candidates in our statistics. 

• NGC23 62 (Figure [T51l 

The 5 Myr (|Mavne fc Navlod l2008l ) old cluster NGC2362 
has few spectroscopically-classified objects and/ or few 
detected sources in the 24 /xm band. 60.0% of the sources 
are still primordial optically-thick while none are classified 
as primordial ultra-settled discs or have lost their discs 
already. 30.0 % are clearing their discs via an IOC process. 

• OBlhf (Figure [T6)) 

61.5 % of the discs in the 5 Myr (|Hernandez et al.l 120061 ) 
old star-forming region OBlbf are still primordial optically- 
thick discs and one disc (7.7%) is classified as primordial 
ultra-settled. 30.8 % of the sources are clearing via the IOC 
process but no sources have lost their discs so far. 

• Upper Sc o (Figure p?]! 

The 5 Myr (|Hernandez et al.l l2007bl ) old Upper Sco re- 
gion seems to be more evolved than most of the other 
star-forming regions considered in this work. Only 9.0 % 
and 2.3 % of the discs, respectively, are still primordial 
optically-thick or classified as primordial ultra-settled. 
80.5 % have already lost their discs and 6.8 % are clearing 
their discs via the IOC process. Only one disc (0.8 %) is 
found at the edge of the HD regime. 

• rj Cha (Figure [T51> 

For the 6 Myr (|Mamaiek et al.l Il999l ) old star-forming 
region rj Cha 35.7% of the objects have been classified as 



primordial optically-thick. 57.1% have already depleted 
their discs, while one disc (7.1 %) is undergoing IOC. 

• 25 On (Figure [T9ll 

In the 7 - 10 Myr (|Briceno et al.ll2007l ) old star-forming re- 
gion 25 Ori objects have barely been detected in the 24 jim 
band. 37.5 % of the objects have primordial optically-thick 
discs, whereas one source (12.5%) has cleared its disc al- 
ready. 37.5 % of the objects are clearing their discs from the 
inside-out. 

3.2 Time-scales 

With the sample of YSOs in different star-forming regions 
becoming larger and larger, it is now possible to estimate the 
typical time-scales for the disc dispersal phase, even though 
cluster ages of course always introduce a large uncertainty 
in the estimates. In order to estimate the time-scale of the 
disc dispersal phase on the basis of our sample, we have 
considered star-forming regions that have more than 10 ob- 
jects. We combine boxes A and C into a single 'primordial' 
category and boxes D and E into a single 'evolved discs' 
category. All objects in box B remain classified as disc-less 
sources. We count objects to the left of the A regime in the 
'primordial' category, if K- [8]>1.5. 

The disc-evolution time-scale of a star-forming region 
can be estimated by r = age-^^*^^^^, where age is the esti- 
mated age of the star-forming region and Ntot and Nevoived 
are the numbers of the total and evolved disc populations, 
respectively. In Table|6]we list the individual time-scales for 
the populations of four spectral type intervals T{SpT) as well 
as their average (t{SpT)) and the cluster time-scale Tc. The 
absolute errors arise from the error propagation of the stan- 
dard statistical errors (e.g. y/N) and the uncertainty in age. 
Note that we arbitrarily assign an error of 10 % to those 
cases where no error is given in the literature. 

Typical transition time-scales for the considered star- 
forming regions are of order 10^ yrs. IC348, OBlbf, 
NGC2264 and Lup HI seem to show larger transition time- 
scales of approximately 1.3 Myrs. The average time-scale 
across all spectral type {(r(SpT))) = 6.9 ■ 10"" yrs is roughly 
the same as the average cluster time-scale (tc) = 6.6-10^ yrs. 
This is partially because there is no significant difference for 
r(SpT) amongst spectral types. 

We further illustrate that disc dispersal time-scales ap- 
pear to be independent of spectral type. We plot in Fig- 
ure 1201 the time-scale ratio r(ii)/r(iii) (K and M stars) 
and show that this ratio is consistent with unity. This sug- 
gests that there is no significant dependence of the time- 
scale on stellar mass, as has already been pointed out by 
lErcolan~ et all ((201 la! ) , who performed a spatial analysis of 
the distribution of K and M stars with discs in young star- 
forming regions and found no significant difference in the 
distributions. 

In Figure 1211 we show histograms, correlating the age 
of the star-forming regions with the number of objects in 
each evolutionary stage. The expected correlation of per- 
centage of evolved and disc-less objects with age is not clear 
from the plot. As well as the usual uncertainty in the ages, 
a number of other factors may contribute to the scatter. 
Very importantly, identification of class HI sources, crucial 
for the statistics, relies on X-ray observations, to make unbi- 
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ased classifications. Unfortunately, only 4 out of 15 regions 
have been sufficiently backed up by X-ray observation. 7 of 
the regions contain some classified X-ray sources, while the 
remaining 4 used different criteria for their member classifi- 
cation. 



4 SUMMARY & CONCLUSIONS 

We have calculated the SEDs of protoplanetary discs of dif- 
ferent spectral types, geometries, settling and inclination. 
We then considered the evolution of the infrared colours (K - 
[8] vs. K-[24]) of the model disks as they disperse accord- 
ing to different scenarios (homologous depletion, inside-out 
and outside-in clearing). Based on our models we propose 
a new diagnostic infrared colour-colour diagram to classify 
the evolutionary stage of YSOs. We have applied our in- 
frared colour-colour diagnostic diagram to classify YSOs in 
15 nearby star-forming regions and study the evolution of 
their disc populations. We estimate time-scales for transi- 
tion phase of typically a few 10^ years independent of stellar 
mass. 

We conclude that, regardless of spectral type, current 
observations point to fast inside-out clearing as the preferred 
mode of disc dispersal. 
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Figure 2. Evolution of a razor-thin disc via HD, IOC and OIC for all four spectral type intervals as described in Section [2. II and analysed 
in Section 12.3.11 Point F are face-on (cosi = 0.95) and point E are edge-on (cosi = 0.05) disc inclinations. The points between the F 
and E represent interme diate disc inclinations in steps of A cosi = 0.10. The red point P represents the approximate locus of the stellar 
photosphere as given bv lLuhman et ah I ||2010|) . The stars mark the HD tracks for face-on and edge-on inclinations (green and blue lines, 
respectively). The triangles mark the IOC tracks for face-on and edge-on inclinations (yellow and cyan lines, respectively). The diamonds 
mark the OIC tracks for face-on and odgc-on inclinations (black and purple lines, respectively). 
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Figure 3. K-[8] vs. K-[24] evolutionary diagnostic diagrams (dashed lines represent objects with envelopes). 
The values of the boundaries are given in the Appendix. 
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Figure 4. Disc evolution diagnostic diagram (dashed lines represent objects with envelopes) applied to the YSOs in the 1 Myr old 
cluster Taurus. (A: 51.2% primordial optically-thick, C: 0.0% primordial ultra-settled, B; 23.2% disc-less, D: 13.7% insidc-out clearing, 
E: 0.5 % homogeneous draining) 20 sample points lie outside the limits of this plot, but they are still included in the final statistics. 
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Figure 5. Disc evolution diagnostic diagram (dashed lines represent objects with envelopes) applied to the YSOs in the 1 Myr old cluster 
Taurus. (A: 41.6% primordial optically-thick, C: 0.8% primordial ultra-settled, B: 31.5% disc-less, D: (17.9) 19.1% insido-out clearing, 
E: (2.3) 3.5% homogeneous draining) 5 sample points lie outside the limits of this plot, but they arc still included in the final statistics. 
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Figure 6. Disc evolution diagnostic diagram (dashed lines represent objects with envelopes) applied to the YSOs in the 2 Myr old cluster 
NGC2068/71. (A: 81.5% primordial optically-thick, C: 0.0% primordial ultra-settled, B: 0.0% disc-less, D: 14.8% inside-out clearing, 
E: 0.0% homogeneous draining) 
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Figure 7. Disc evolution diagnostic diagram (dashed lines represent objects with envelopes) applied to the YSOs in the 2 Myr old 
cluster NGC2264. (A: 25.2% primordial optically-thick, C: 0.6% primordial ultra-settled, B: 9.9% disc-less, D: (56.2) 57.2% inside-out 
clearing, E: (0.3) 1.3% homogeneous draining) 
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Figure 8. Disc evolution diagnostic diagram (dashed lines represent objects with envelopes) applied to the YSOs in the 2-3 Myr old 
cluster Cha I. (A; 75.9% primordial optically-thick, C; 0.0% primordial ultra-settled, B: 0.0% disc-less, D: 3.4% inside-out clearing, 
E: 3.4% homogeneous draining) 1 sample point lies outside the limits of this plot, but it is still included in the final statistics. 
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Figure 9. Disc evolution diagnostic diagram (dashed lines represent objects with envelopes) applied to the YSOs in the 2-3 Myr old 
cluster IC348. (A; 23.8% primordial optically-thick, C: 0.8% primordial ultra-settled, B: 3.8% disc-less, D: (69.6) 70.0% inside-out 
clearing, E: (0.0) 0.4% homogeneous draining) 3 sample points lie outside the limits of this plot, but they arc still included in the final 
statistics. 
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Figure 10. Disc evolution diagnostic diagram (dashed lines represent objects with envelopes) applied to the YSOs in the less than 3 
Myr old cluster NGC1333. (A: 66.7 % primordial optically-thick, C: 0.0 % primordial ultra-settled, B: 0.0 % disc-lcss, D: 28.2 % insidc-out 
clearing, E: 0.0% homogeneous draining) 9 sample points lie outside the limits of this plot, but they arc still included in the final 
statistics. 
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Figure 11. Disc evolution diagnostic diagram (dashed lines represent objects with envelopes) applied to the YSOs in the 4 Myr old 
cluster Tr 37. (A: 89.5% primordial optically-thick, C; 0.0% primordial ultra-settled, B: 0.0% disc-less, D; 7.0% inside-out clearing, 
E: 0.0% homogeneous draining) 
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Figure 12. Disc evolution diagnostic diagram (dashed lines represent objects with envelopes) applied to the YSOs in the 2-6 Myr old 
cluster Serpens. (A: 40.6% primordial optically-thick, C: 1.4% primordial ultra-settled, B: 18.8% disc-less, D: (10.1)11.6% inside-out 
clearing, E: (10.1) 11.6 % homogeneous draining) 2 sample points lie outside the limits of this plot, but they arc still included in the final 
statistics. 
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Figure 13. Disc evolution diagnostic diagram (dashed lines represent objects with envelopes) applied to the YSOs in the 2-6 Myr old 
cluster Lup III. (A: 39.0% primordial optically-thick, C: 1.2% primordial ultra-settled, B: 31.7% disc-less, D: 14.6% inside-out clearing, 
E: 7.3 % homogeneous draining) 2 sample points lie outside the limits of this plot, but they are still included in the final statistics. 
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Figure 14. Disc evolution diagnostic diagram (dashed lines represent objects with envelopes) applied to the YSOs in the 4-5 Myr old 
cluster Cha II. (A: 63.0% primordial optically-thick, C: 0.0% primordial ultra-settled, B: 25.9% disc-less, D: 7.4% insidc-out clearing, 
E: 0.0% homogeneous draining) 1 sample point lies outside the limits of this plot, but it is still included in the final statistics. 
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Figure 15. Disc evolution diagnostic diagram (dashed lines represent objects with envelopes) applied to the YSOs in the 5 Myr old 
cluster NGC2363. (A: 60.0 % primordial optically-thick, C: 0.0 % primordial ultra-settled, B: 0.0 % disc-less, D: 30.0 % inside-out clearing, 
E: 0.0% homogeneous draining) 
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Figure 16. Disc evolution diagnostic diagram (dashed lines represent objects with envelopes) applied to the YSOs in the 5 Myr old 
cluster OBlbf. (A: 61.5% primordial optically-thick, C; 7.7% primordial ultra-settled, B: 0.0% disc-less, D: 30.8% inside-out clearing, 
E: 0.0% homogeneous draining) 
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Figure 17. Disc evolution diagnostic diagram (dashed lines represent objects with envelopes) applied to the YSOs in the 5 Myr old 
cluster Upper Sco. (A: 9.0% primordial optically-thick, C: 2.3% primordial ultra-settled, B; 80.5% disc-less, D: 6.8% inside-out clearing, 
E: 0.8 % homogeneous draining) 
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Figure 18. Disc evolution diagnostic diagram (dashed lines represent objects with envelopes) applied to the YSOs in the 6 Myr old 
cluster r] Cha. (A: 35.7% primordial optically-thick, C: 0.0% primordial ultra-settled, B: 57.1% disc-less, D: 7.1% inside-out clearing, 
E: 0.0% homogeneous draining) 
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Figure 19. Disc evolution diagnostic diagram (dashed lines represent objects with envelopes) applied to the YSOs in the 7-10 Myr old 
cluster 25 Ori. (A: 37.5% primordial optically-thick, C: 0.0% primordial ultra-settled, B: 12.5% disc-less, D: 37.5% inside-out clearing, 
E: 0.0% homogeneous draining) 
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Taurus Taurus NGC2264 Cha I IC348 NGC1333 Serpens Lup III Cha II 0B1bf Upper Sco 25 Ori 
(Luhman) (Rebull) 

Star-forming regions 

Figure 20. Timc-scalc ratio t{II)/t{III) for K and M stars for the star-forming regions considered in this work with Ntot > 10. Tr37, 
NGC2068/71 and rj Cha arc not listed, because they lack evolving objects in one or both spectral type intervals. 
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Figure 21. Number of objects in each evolutionary stage (as explained in the legend) as a function of the cluster age. 
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APPENDIX A: 



We list here the boundary points of the evolutionary boxes 
to allow automatized application to observational data. 



The boundaries for primordial discs with envelopes. 
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II 
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II 
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II 
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[2.99, 3.14, 3.64, 3.86, 4.48, 4.43, 3.99, 3.78, 3.14, 2.99] 

[3.38, 3.28, 2.34, 1.37, 0.76, 1.03, 1.84, 2.85, 1.84, 3.14] 

[2.6,2.33,1.53,1.16,1.02,1.34, 3.] 

[0.77, 0.87, 1.1, 1.16, 1.2, 1.12, 0.91, 0.86, 0.73, 0.77] 

[1.53, 1.53, 1.08, 0.62, 0.37, 0.29, 0.54, 1.02, 0.54, 0.73] 

[7., 6.77, 5.61, 4.18, 2.95, 2.95, 3.12, 5.3] 

[2.71, 3.22, 3.87, 4.17, 4.57, 4.18, 3.7, 2.83, 2.71] 

[3.44, 2.69, 1.84, 0.98, 1.31, 2.05, 2.95, 2.05, 2.56, 3.22] 

[2.8, 2.6, 1.67, 1.33, 1.17, 1.19, 1.41, 3., ] 

[0.83, 0.82, 0.91, 1.01, 1.28, 1.33, 1.16, 0.89, 0.83] 

[1.65, 1.34, 0.85, 0.47, 0.37, 0.68, 1.19, 0.68, 0.69, 0.82] 



ENV^ ' 
ENVf^ ' 


[24] ^ 
-[8] ^ 


[3.14, 5.75,6.8,8.9; 
[0.73, 0.91,1.28,3.; 


ENV^j ' 
ENVl^j ' 


[24] ^ 
-[8] ^ 


[3.22,5.3,6.9,9.] 
[0.82,0.85,1.38,3.; 


ENV^,- 
ENV;",, 


[24] ^ 
-[8] ^ 


[3.73,4.75,6.7,8.9; 
[0.94,1.0,1.65,3.] 


ENVry' 

ENVl'y 


[24] ^ 
-[8] ^ 


[3.95, 5.2, 7., 8.5] 
[1.23, 1.27,2.05,3.; 



^iii 


[24] 


= [7., 6.2, 5.01, 2.96, 3.01, 6.] 


'-^iii 


[24] 


= [2.96, 2.97, 3.69, 3.87, 4.3, 3.73, 2.96] 


^UI 


[24] 


= [3.45, 3.33, 2.66, 1.33, 1.63, 2.16, 2.58, 2.96] 


^iii 


-[8] 


= [2.6,2.1,1.41,1.11,1.3,3.] 


^III 


-[8] 


= [0.91, 0.97, 1.28, 1.33, 1.32, 0.94, 0.91] 


^ni 


-[8] 


= [1.55, 1.56, 1.31, 0.67, 0.55, 0.81, 0.83, 0.91] 



Afv' = [7., 6.46, 5.4, 4.11, 3.22, 3.2, 3.47, 5.8] 

^K- [24] ^ J2.97, 3.95, 4.42, 3.88, 3.08, 2.97] 

E^y'^'^^^ = [3.67,3.39,2.6,2.12,1.48,1.6,2.45,2.76,3.02,3.2] 

Afy~^^^ = [2.8,2.43,1.74,1.51,1.35,1.37,1.62,3.] 

C^"'*' = [1.13,1.23,1.59,1.54,1.21,1.13] 

= [1.75,1.74,1.29,1.05,0.77,0.72,1.04,1.06,1.08,1.15] 



From the IOC runs for finite-thickness discs, we can 
conclude further that all objects found in the D box with 
colours bluer in K - [24] than 4 can be considered as evolved 
ultra-settled objects. Objects redder in K- [24] than 5 
represent the evolved discs which are either settled or 
mixed and/ or flared. 



